Non-Hermitian non-equipartition theory for trapped particles

The equipartition theorem is an elegant cornerstone theory of thermal and statistical physics. However, it fails to address some contemporary problems, such as those associated with optical and acoustic trapping, due to the non-Hermitian nature of the external wave-induced force. We use stochastic calculus to solve the Langevin equation and thereby analytically generalize the equipartition theorem to a theory that we denote the non-Hermitian non-equipartition theory. We use the non-Hermitian non-equipartition theory to calculate the relevant statistics, which reveal that the averaged kinetic and potential energies are no longer equal to kBT/2 and are not equipartitioned. As examples, we apply non-Hermitian non-equipartition theory to derive the connection between the non-Hermitian trapping force and particle statistics, whereby measurement of the latter can determine the former. Furthermore, we apply a non-Hermitian force to convert a saddle potential into a stable potential, leading to a different type of stable state.

The equipartition theorem is an elegant cornerstone theory of thermal and statistical physics.However, it fails to address some contemporary problems, such as those associated with optical and acoustic trapping, due to the non-Hermitian nature of the external wave-induced force.We use stochastic calculus to solve the Langevin equation and thereby analytically generalize the equipartition theorem to a theory that we denote the non-Hermitian nonequipartition theory.We use the non-Hermitian non-equipartition theory to calculate the relevant statistics, which reveal that the averaged kinetic and potential energies are no longer equal to k B T/2 and are not equipartitioned.As examples, we apply non-Hermitian non-equipartition theory to derive the connection between the non-Hermitian trapping force and particle statistics, whereby measurement of the latter can determine the former.Furthermore, we apply a non-Hermitian force to convert a saddle potential into a stable potential, leading to a different type of stable state.
Brownian motion is a fundamental type of thermal motion that is of paramount importance in various scientific and technological applications.The underlying mechanism of Brownian motion was explained by Einstein in 1905 1 and Smoluchowski in 1906 2 , and its characteristic jittery movement results from the irregular bombardment of fluid molecules.This movement is random, complex, and unrepeatable, and thus Brownian motion can only be described by statistical theoretical treatments.Perhaps the most important and elegant result of studies on Brownian motion is the equipartition theorem (ET), which was first devised in 1843 for an equipartition of kinetic energies 3 , and subsequently generalized 4 , such that it became a cornerstone of classical statistical physics 5 .The ET states that at thermal equilibrium, every quadratic degree of freedom has an average energy of k B T/2, where k B is the Boltzmann constant and T is temperature.Furthermore, the ET plays a significant role in describing a broad range of physical scenarios, including the ideal gas law 5 , the Dulong-Petit law for specific heat capacities of solids 6 , Graham's law of effusion 7 , the extreme relativistic ideal gas in astrophysics 5,8 , and Johnson-Nyquist noise 9 .
However, the ET is inadequate for certain contemporary problems, such as those associated with optical/acoustic trapping [10][11][12][13][14][15][16][17][18] and binding [19][20][21][22][23][24][25][26][27][28][29][30][31][32] .These involve a single particle or a system of particles bound at mechanical equilibrium by an external optical or acoustic wave, which causes the particle or system to be out of thermal equilibrium.The wave exchanges energy with the particle(s), leading to the generation of nonconservative forces that drive the particle or system into a non-Hermitian state 31 .This drastically alters the physical characteristics of a particle or system, and thus the ET fails to describe a particle or system whenever non-Hermitian manipulation is performed.External nonconservative forces compete with ambient damping, with the former pumping energy into particles and the latter removing energy from particles.Regarding optical manipulation, non-Hermitian forces are expected to play a vital role in optical trapping under vacuum (low damping) and a direct and significant role in optical trapping in air under atmospheric conditions (intermediate damping).Moreover, during optical trapping in water (heavy damping), non-Hermitian forces are expected to drive particles' characteristic vibrational modes into a non-orthogonal state that deviates from Hermitian physics.
We use stochastic calculus to solve the Langevin stochastic differential equation with nonconservative trapping forces and thereby generalize the ET to the non-Hermitian non-equipartition (NHNE) theory (may be equivalently termed non-Hermitian non-equipartition theorem).We also use the Verlet algorithm 63 to validate the analytical stochastic calculations through numerical simulations.Our results significantly deviate from the k B T/2 predicted by the ET, and the energies are no longer equipartitioned among the different degrees of freedom.Moreover, we observe qualitative discrepancies.For instance, a sufficiently large nonconservative force can destabilize an otherwise stable trap 14,31 , whereas a particle in an originally unstable saddle potential can be stabilized by an appropriate nonconservative force.The NHNE theory is also capable of measuring forces other than conservative trapping forces 64 , as it can measure non-Hermitian forces and the repulsive forces of a particle in a saddle potential.Additionally, we discuss the NHNE theory for N > 1 Brownian particles, which qualitatively captures accurate numerical results by applying an approximate analytical theory.

Statement of the NHNE theory
Consider a spherical particle of mass m that is immersed in a fluid and is confined within the vicinity of an equilibrium r 0 = ðx 0 ,y 0 ,z 0 Þ by a force field FðrÞ = ðF x ,F y ,F z Þ.The dynamics of the particle are governed by the Langevin stochastic differential equation 65 : where r = ðx,y,zÞ represents the position of the particle; γ = 6πηa is the friction coefficient; a is the particle radius; η is viscosity; and AðtÞ = ðA x ðtÞ,A y ðtÞ,A z ðtÞÞ, which obeys the fluctuation-dissipation theorem, is the Gaussian-distributed random force due to Brownian fluctuations.AðtÞ has correlations of hA i ðtÞA j ðt 0 Þi = 2γk B Tδ i,j δðt À t 0 Þ, where A i ðtÞ represents the ith component of AðtÞ, δ i,j is the Kronecker delta function, and δðtÞ is the Dirac delta function.
Near the equilibrium position r', where Fðr 0 Þ = 0, FðrÞ in Eq. ( 1) may be sufficiently approximated by its linear term (as the zeroth-order term vanishes), as follows: where calculated at r = r 0 is the force matrix 22,31 , where F i and r i are the ith components of FðrÞ and r, respectively.$ K ij and F i can be calculated numerically.For optical trapping, which we consider later, the force calculation is performed using Mie scattering theory and the Maxwell stress tensor 22,31,66 .It is equally valid to define the force matrix with an extra minus sign as FðrÞ ≈ À $ K Á ðr À r 0 Þ, where a positive eigenvalue of $ K indicates a stable mode.However, including an extra minus sign is a matter of preference and has no physical significance.Without loss of generality, one could take the equilibrium to be the origin, i.e., r 0 = 0.By solving Eqs. ( 1) and (2) using stochastic calculus (Supplementary Notes 1, 2), we obtain where r i (v i ) denotes the ith component of the displacement (velocity), the columns of  3) is one of the main results of this paper and generalizes the ET to non-Hermitian systems.A detailed derivation is presented in Supplementary Notes 1, 2. In general, unlike in the ET, the average energies in a non-Hermitian system associated with different degrees of freedom are not equipartitioned by Eq. (3).Hence, the generalization of the ET theory that it represents is referred to as the NHNE theory.We note that the NHNE theory, denoted by Eq. ( 3), can be applied to any force matrix $ K, including both Hermitian and non-Hermitian matrices, for a trapped Brownian particle.
The NHNE theory for optical trapping We provide a concrete example by considering optical trapping as a non-Hermitian system.In the case of a particle trapped by a wave, the non-Hermiticity of the force matrix results from light scattering 31,57,58 .We utilize Eq. ( 3) to investigate the specific scenario of a single particle trapped by light; the NHNE theory is also applicable to other mechanical systems, including acoustic trapping and binding systems.The motion of the trapped particle along the z-axis of a typical trapping beam is independent of the transverse motions whenever 31 The conditions sufficient for Eq. ( 4) to hold include but are not limited to (i) the equilibrium state of the system exhibits mirror symmetry about the z = z 0 plane (e.g., in the case of identical counterpropagating beams), (ii) the system exhibits mirror symmetries about the x = x 0 and y = y 0 planes at the equilibrium position (e.g., in the case of a linearly polarized Gaussian beam), or (iii) the system exhibits rotational symmetry (e.g., with a circularly polarized Gaussian beam).Here, ðx 0 ,y 0 ,z 0 Þ is the equilibrium position of the particle.Equation ( 4) enables us to resolve the transverse motion by using the reduced twodimensional (2D) non-Hermitian force matrix , where By implementing a coordinate transformation that diagonalizes the symmetric part of $ K 0 2D , we obtain, without loss of generality, , where $ R is a rotation matrix.Here, the existence of the anti-symmetric component g in the force matrix is a consequence of non-Hermiticity.One can infer from the force field produced by g that this nonconservative force revolves the particle around the mechanical equilibrium.
Then, Eq. ( 3) is simplified to NHNE theory, which is fundamentally distinct from the ET.Here, we note that the temperature T refers to the ambient temperature.It is unchanged for transparent particles since the work done by the non-Hermitian force is on the order of k B T (comparable with the energy of a single molecule in the ambient), thus completely negligible for the entire ambient.When g = 0, $ K 2D is Hermitian, and Eq. ( 5) reduces to the ET, where 1 2 mhv 2 However, when g ≠ 0, Eq. ( 5) can significantly deviate from the ET.The equality of kinetic energy in Eq. ( 5) results from our aiming to simplify the equation by selecting a specific coordinate system (where ).The general form for arbitrary Generally, the kinetic energies along two orthogonal transverse directions are not equal, as explained in Supplementary Note 3. If the non-Hermitian force matrix is scaled as , the averaged quantities in Eq. ( 5) are nonlinear with respect to the scalar P, as altering P changes the ratio between the trapping forces and the random forces.This also applies to Eq. (3).In optical and acoustic trapping, P is proportional to the incident power, as detailed in Supplementary Note 6.
We now apply Eq. ( 5) to investigate a Brownian particle (with refractive index n = 1.57and radius a = 0.5 μm) that is illuminated and trapped by two counter-propagating Gaussian beams that each have an input power of 1.0 mW.Such a system exhibits a non-Hermitian force matrix that is tunable by varying the incident polarization p = x cosðζ Þ + iŷ sinðζ Þ.A similar system was examined in ref.31, but without considering Brownian motion.The values of the components of $ K 2D versus the incident polarizations are illustrated in Fig. 1a, while the corresponding eigenvalues, are shown in Fig. 1b, where k xx and k yy represent the restoring force constants in the x and y directions, respectively.The parameter g represents the torque driven by the orbital angular momentum of the elliptically or circularly polarized Gaussian beam, which is created by the spin-orbit angular momentum conversion during focusing by the objective lens 67 .As ζ is continuously tuned from 0°to 45°, the polarization gradually varies from linear (Fig. 1c) to elliptical and then to circular (Fig. 1d).With linear polarization (ζ = 0°), g = 0 due to reflection symmetries on the xz and yz planes, and thus the system is effectively Hermitian.In addition, k xx ≠ k yy due to polarization aberrations 68 .As ζ increases to (for example) 10°, the left circular polarization becomes stronger than the right circular polarization, which breaks the balance of the orbital angular momentum originally The focused beam has a numerical aperture of 0.9 and a filling factor of 1.0, and the beam power for each beam is normalized to 1.0 mW.For different viscosities, η = 3.0 μPa•s (e-g) and η = 18.4 μPa•s (h-j), average energies 1 2 mhv 2 x i and 1 2 mhv 2 y i, À 1 2 k xx hx 2 i and À 1 2 k yy hy 2 i, and À 1 2 ghxyi are presented as lines (Anal.: analytical values) and circles (Num.: numerical values).
presented in the linear polarization and generates a finite g.Moreover, the focused spot approaches cylindrically symmetric, and thus the difference between k xx and k yy decreases.At ζ = 21.6°, which is indicated by the red dashed line in Fig. 1b, an exceptional point (EP) emerges, where the real parts of the eigenvalues merge and their imaginary parts split.The EP arises due to a switch in the sign of the quantity under the square root in Eq. ( 6), which causes the originally real eigenvalues to become complex.This is expected, as ðk xx À k yy Þ 2 ultimately reaches zero (at ζ = 45°) due to rotational symmetry, while g increases monotonically from zero.As a result, the existence of an EP at g 2 = ðk xx À k yy Þ 2 is inevitable, regardless of the specifics of a system.
The values of K ± in Eq. ( 6) are real and negative (indicating stability) to the left of the EP, whereas they become complex (indicating instability if ambient damping is insufficient) 31 to the right of the EP. Figure 1e shows the kinetic energies of 1 2 mhv 2 x i and 1 2 mhv 2 y i plotted for a viscosity of η = 3.0 μPa•s, whereas Fig. 1f shows the quantities À 1 2 k xx hx 2 i and À 1 2 k yy h y 2 i.Crucially, the last two quantities, referred to as potential energies, only represent the conservative interaction, but not the entire non-Hermitian interactions.The latter cannot be expressed as the gradient of a potential.According to the ET, all four quantities should equal k B T/2.However, they deviate from k B T/2, except when the system is effectively Hermitian under linear polarization (ζ = 0°).Figure 1g also displays À 1 2 ghxyi, which has non-zero values, indicating the presence of non-Hermitian couplings between the x and y motions.Due to the axial symmetry, À 1 2 ghxyi = 0 under circular polarization (ζ = 45°), despite non-Hermitian coupling still being present.We compare the results obtained using the NHNE theory (Eq.( 5)) with those generated by a numerical Verlet simulation (marked by circles in Fig. 1e-g).The Verlet simulation utilizes the "exact" optical force field computed from Mie scattering theory 66 , instead of the linearized force based on the force matrix.In the Mie scattering theory, the expansion series are truncated at L max = ka + 4ðkaÞ 1=3 + 2 69,70 , with k = 2π λ being the wavenumber.In some cases, we verified the convergence of our calculations by comparing the calculation truncated at L max with those truncated at L max + 5. Further details on the Verlet simulation can be found in the Methods.The results show remarkable agreement, with small deviations found only near areas of circular polarization.This is because circular polarization has the largest g, which leads to the fastest rate of energy pumping, allowing a particle to move further away from the origin than under other types of polarization, such that the linear approximation in Eq. ( 2) becomes less accurate.The agreement with the "exact" force calculation demonstrates the sufficiency of the linear approximation used to derive Eq. ( 5).
Figure 1h-j show the same settings as Fig. 1e-g, except that the former has a viscosity η = 18.6 μPa•s (air), whereas the latter has η = 3.0 μPa•s.As the damping increases, the maximum deviations from the ET decrease from ~200% to ~20%.This indicates that although the non-Hermiticity is significantly suppressed, it cannot be ignored, even for optical trapping in air.A large dissipation implies that the energy the particle receives from light is dissipated quickly, preventing its accumulation.Regardless of the level of damping, the presence of non-Hermitian forces ensures that the vibrational eigenmodes are always non-orthogonal, indicating that the non-Hermiticity of the force matrix cannot be ignored in any case.In all cases, the results from the NHNE theory remain accurate.The non-Hermiticity of optical trapping also depends strongly on the particle radius (a), as detailed in Supplementary Note 8.
Interestingly, À 1 2 k xx hx 2 i and À 1 2 k yy h y 2 i in Fig. 1f, i are not the same.In fact, in Fig. 1i, À 1 2 k yy h y 2 i can even be less than k B T/2.This observation can be explained by the virial theorem 71 , 1 2 mhv 2 i i = K ij hr i r j i, which remains applicable even in the non-Hermitian case (Supplementary Note 3).In addition, according to the virial theorem, one can obtain 1 2 mhv 2 Figure 2 shows the average rate at which work is done on the particle by the dissipative force Àhγv Á vi, random force hAðtÞ Á vi, and optical force hFðrÞ Á vi, calculated analytically (lines) and numerically (circles), as a function of polarizations (Fig. 2a) and viscosities (Fig. 2b).As required by energy conservation, we have Àhγv Á vi + hAðtÞ Á vi + hFðrÞ Á vi = 0, as detailed in Supplementary Note 4. As AðtÞ is a random force that does not depend on v, its power is expected to be the same for a particle moving freely and a particle moving under a force field, resulting in a constant hAðtÞ Á vi = 2γk B T m in Fig. 2a.Moreover, in Fig. 2a, hFðrÞ Á vi increases with the non-Hermiticity of the force matrix as ζ increases, which accounts for the increased energies in the NHNE theory results shown in Fig. 1e-j.As η increases, the role of the optical force initially increases and then decreases to zero, because at a large η, work is done at a high rate by the Brownian fluctuation hAðtÞ Á vi and the corresponding motion is heavily damped.Surprisingly, the optical force does no work in a perfect vacuum, but it does work when fluctuation and damping are present; this highlights an unexpected role of Brownian motion in light-driven machines.
The NHNE theory for a saddle potential with a nonconservative force field A non-Hermitian force matrix consists of a potential energy term (the symmetric part of $ K 2D ) and a nonconservative term (the antisymmetric part of $ K 2D ).A saddle potential traps a particle in one direction (when k xx < 0) and repels the particle in the other direction (when k yy > 0), making it unstable.Counter-intuitively, a non-Hermitian force can stabilize a particle in a saddle potential if the trace of 3c) to circular (ζ = 45°, Fig. 3d), the matrix elements of $ K 2D , namely k xx , k yy , g, and -g, are illustrated with lines in Fig. 3a for optical trapping of a hollow sphere (with refractive index n = 1.57, inner radius 0.49 μm, and outer radius a = 0.70 μm), where the saddle potential (defined by k xx Á k yy < 0) occurs on the left-hand side of the black dashed line.Figure 3b plots ReðK ± Þ and ImðK ± Þ versus ζ, wherein the EP is marked with a red dashed line, and the neutral point (NP) with K + = 0 and K À <0 is marked with a blue dashed line.Coincidentally, the black dashed line (Fig. 3a) and the EP are very close.On the left-hand side of the NP, the positive K + mode repels the particle, regardless of the ambient damping level, whereas on the right-hand side of the EP, the complex modes K ± destabilize the particle when the background damping is not sufficiently large.Unexpectedly, between the NP and the EP, the particle is always stable.Similar to Fig. 1, we depict the NHNE in Fig. 3e-j) with viscosity η = 3.0 μPa•s (η = 18.4 μPa•s) versus polarizations.At the NP, À 1 2 k xx hx 2 i, À 1 2 k yy h y 2 i, and À 1 2 ghxyi diverge because there is no confinement by optical force along the neutral direction.However, the kinetic energies remain finite.For η = 3.0 μPa•s, hv 2 x i, hv 2 y i, hx 2 i, h y 2 i, and hxyi diverge beyond the polarization characterized by ζ = 41.7°,i.e., after the EP, because K ± are complex and the background damping is insufficient, making the trapping unstable.Figure 3e-j present the analytical (lines) and numerical (circles) calculations.Inconsistencies are only observed when the particle is far away from the equilibrium position.As k yy varies from positive to negative due to the increase in ζ, À 1 2 k yy h y 2 i also varies from negative to positive (Fig. 3f, i), i.e., from repulsion to trapping.This is a unique phenomenon for a saddle potential with a non-Hermitian force field.As an additional example of saddle potential, we also consider a uniform layer of dielectric coated on a gold sphere for the purpose of optical trapping, as illustrated in Supplementary Note 9.

Trajectories for optical trapping
The trajectories of particles trapped by non-Hermitian forces and subject to Brownian motion are depicted in Fig. 4a, b (trapping potential) and Fig. 4c, d (saddle potential), with their initial positions marked by black dots.For the trapping potential at η = 1.0 μPa•s (Fig. 4a), the viscosity is too low to confine the particle after the EP.However, when the viscosity is increased to η = 3.0 μPa•s, as shown in Fig. 4b, the Brownian particle is stable at all of the polarizations.We can still observe the expansion of the trajectories as ζ increases 58 , which is a result of the increasing hx 2 i and h y 2 i shown in Fig. 1f.In addition, the trapped particle exhibits Brownian fluctuations near the equilibrium, and its orientation varies with ζ (see Supplementary Note 5).
For the saddle potential, the particle escapes before the NP (ζ < 26.9°), due to the repulsive force (Fig. 4c, d), and after ζ = 32.6°(ζ = 41.7°), when the viscosity is as low as η = 1.0 μPa•s (η = 3.0 μPa•s), as shown in Fig. 4c (Fig. 4d).On the NP, the Brownian particle is trapped in only one direction, whereas the fluctuation force repels the particle far away in another direction.Between the NP and the EP (26.9°< ζ < 32.0°), the Brownian particle is trapped stably, even though the optical force repels the particle in one direction.The videos for the three-dimensional trajectories for each case depicted in Fig. 4a-d are available in Supplementary Movies 1-4, and detailed phase diagrams for both trapping and saddle potentials can be found in Supplementary Note 6.

The NHNE theory for multiple particles
Compared with the NHNE theory for a single particle, the NHNE theory for multiple particles is significantly more complex due to hydrodynamic interactions between particles 72 .Currently, there is no exact Furthermore, an approximate solution can be derived if the hydrodynamic interactions are neglected (see Supplementary Note 7): The approximate analytical solutions semi-qualitatively agree with the numerical results.
We consider a linear chain of spheres with N = 2 (Fig. 5a) that are optically trapped and bounded in a low vacuum by two coherent linearly polarized plane waves (η = 1.0 μPa•s) propagating along the chain axis.The two plane waves have different intensities I 1 and I 2 (I 1 = 9I 2 ), respectively, creating an unbalanced propagation that favors the non-Hermitian force fields 73,74 .Each sphere has a radius of a = 0.2 μm and a refractive index of n.We search for an equilibrium configuration for the pair of spheres and consider their motion along the z-axis.The red lines in Fig. 5b (with separation D ≈ λ and refractive index n = 1.1) and Fig. 5c (with D ≈ 4λ and n = 1.2) plot the averaged kinetic energies 1 2 mhv 2 z,i i, which are calculated using Eq. ( 7) for each sphere i versus the averaged intensity I 0 = 1 2 ðI 1 + I 2 Þ.Here, where $ K 0 is the force matrix at I 0 = 1:0W=m 2 .As the componnets of the non-Hermitian force matrix increase with I 0 , the kinetic energies increasingly deviate from k B T/2 and are no longer equipartitioned.We also conduct numerical simulations using Verlet algorithms that take the hydrodynamic interactions into account.The results are depicted as blue dots in Fig. 5b (for a small D) and Fig. 5c (for a large D).The hydrodynamic interactions have a significant impact on the partitioned energies when D is small, i.e., the blue dots deviate from the red line, as shown in Fig. 5b.However, the deviation diminishes when D is large, as the hydrodynamic interaction is weakened, as shown in Fig. 5c.In sum, the approximate multiple particle NHNE theory (Eq.( 7)), from which hydrodynamic interactions are excluded, can semi-qualitatively predict the averaged energies for each degree of freedom, especially when D is large.

Discussion
In this article, stochastic calculus is applied to generalize the ET to deal with non-Hermitian trapping and binding forces.This  generalized theory, denoted by the NHNE theory, enables the calculation of the average energies of a single particle or a group of trapped or bounded particles, even when the force matrix of a system is non-Hermitian.By "generalize", we meant to extend the original ET to address non-Hermitian problems.We note that the NHNE theory reveals the breaking of universality in the original ET by non-Hermiticity, in the sense that the average energies associated with each degree of freedom are no longer equal and depends on the details of the system.This is a development in the study of Brownian motion and has far-reaching implications for a variety of problems associated with modern technology, including those associated with optical/ acoustic trapping and binding, and with other open mechanical systems.
To provide a concrete illustration of the NHNE theory, we focus on optical trapping, which is one example of a non-Hermitian trapping system.The NHNE theory is applied to analyze both the trapping potential (Fig. 1a) and the saddle potential (Fig. 3a).We propose that by experimentally measuring hx 2 i, h y 2 i, hxyi, hv 2 x i, hv 2 y i, and hv x v y i, Eq. ( 5) can be used to determine both the Hermitian (k xx and k yy ) and non-Hermitian (g) force constants.To our knowledge, such an approach is previously limited by the availability of pertinent theories, as a result, they can only be used to perform such indirect measurement of force constants in a heavily damped environment .Thus, the NHNE theory provides a method to directly measure force constants under arbitrary levels of damping, including in vacuum trapping applications.When there is a large damping, some predictions by our theory (such as the average energies) can be very similar to what conventional ET predicts.It might seem like damping is getting rid of non-Hermiticity, but that is not the complete physical picture.For example, the vibrational eigenmodes remain non-orthogonal, irrespective of the damping.Moreover, we make the surprising finding that non-Hermitian forces can stabilize a particle in a saddle potential.Repulsive forces at the microscopic scale can be difficult to measure, due to the absence of a stable equilibrium.We propose that a microparticle located near a saddle potential can be stabilized by non-Hermitian forces.Furthermore, the repulsive force constant in a saddle potential can be determined using the NHNE theory.
We note that the NHNE theory is of relevance to the study of nonreciprocal interaction and active matters [76][77][78][79][80] .These systems are typically complex, involving intricate geometries, and a variety of nonreciprocal and non-Hermitian interactions.Our theory may offer some insights and approximate predictions into these problems.In Eq. ( 3),

Expressions of
where and ω i = ffiffiffiffiffi K i p with K i being the ith eigenvalue of À $ K=m.

Verlet algorithm for Langevin dynamics simulations
We use the Verlet algorithm 63 to solve the stochastic Langevin differential equation in Eq. ( 1) numerically.By definition, the fluctuating force AðtÞ = ðA x ðtÞ,A y ðtÞ,A z ðtÞÞ, which satisfies hA i ðtÞA j ðt 0 Þi = 2γk B Tδ i,j δðt À t 0 Þ, is independent of the particle velocity vðtÞ, and one can assume that AðtÞ is a constant force during each time step h.We denote the fluctuating force during the time interval [t, t + h] ([t, th]) as A + (A À ).The particle positions at time t, th, and t + h are related by rðt + hÞðm + γh=2Þ + rðt À hÞðm À γh=2Þ = 2rðtÞm + h 2 ðFðrðtÞÞ + A + =2 + A À =2Þ, ð10Þ while for the particle velocities, vðt + h=2Þðm + γh=2Þ = vðt À h=2Þðm À γh=2Þ + hðFðrðtÞÞ + A + =2 + A À =2Þ: Here, FðrðtÞÞ denotes the external force (e.g., optical force) exerted on the particles located at rðtÞ.The Verlet results presented in the main text and Supplementary Information are obtained using a time step h = 10 À8 s.The average quantities are based on the average of 10 10 positions or velocities.

Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
in the Methods.Equation (

Fig. 1 |
Fig. 1 | Non-Hermitian non-equipartitioning in optical trapping.a Components of the force matrix $ K 2D = k xx g Àg k yy and b its eigenvalues for an optically trapped dielectric Brownian particle with refractive index n = 1.57and radius a = 0.5 μm in a low vacuum or in air, as the polarization (p = x cos ζ + iŷ sin ζ ) of the incident beams (with wavelength λ = 1.064 μm) varies from linear (ζ = 0°, panel c) to circular

Fig. 2 |
Fig. 2 | Power delivered by different force components.Power delivered by the damping force hγv Á vi (black), the fluctuational force hAðtÞ Á vi (red), and the optical force hFðrÞ Á vi (blue) versus a the polarizations (ζ) and b viscosities (η) in the optical trapping of a Brownian particle (identical to that described in Fig. 1e-g) is presented.The analytical (Anal.)and numerical (Num.)values are denoted by lines and circles, respectively.

Fig. 3 |mhv 2 x i and 1 2
Fig. 3 | Non-Hermitian non-equipartitioning in optical trapping of a saddle potential.Components (a) and eigenvalues (b) of the force matrix $ K 2D = k xx g Àg k yy for an optically trapped dielectric hollow particle with refractive index n = 1.57, inner radius 0.49 μm, and outer radius a = 0.70 μm, in a low vacuum or in air, as the polarization (p = x cos ζ + iŷ sin ζ ) of the incident beams (with wavelength λ = 1.064 μm) varies from linear (ζ = 0°, c) to circular (ζ = 45°, d).The left-hand side of the black dashed line in a corresponds to the saddle potential,

Fig. 5 |
Fig. 5 | Non-Hermitian non-equipartitioning in the optical binding of two spheres.a The two spheres (colored green), each with a radius of a = 0.2 μm and a refractive index of n, are optically bounded by two linearly polarized counterpropagating waves (with wavelength λ = 1.064 μm) in a low vacuum (η = 1.0 μPa•s).The two plane waves have different intensities, with I 1 = 9I 2 , and I 0 = 1 2 ðI 1 + I 2 Þ.The average kinetic energies 1 2 mhv 2 z,i i for the particles (indexed by i, ranging from 1 to 2) versus I 0 are presented as red lines (analytical (Anal.)results excluding hydrodynamic interactions) and blue dots (numerical (Num.)results including hydrodynamic interactions) for different separations (D) and different values of n, namely b D≈λ and n = 1.1, and c D≈4λ and n = 1.2.This analysis focuses on motion along the z-axis.

Fig. 4 |
Fig. 4 | Trajectories for an optically manipulated Brownian particle.Numerically calculated trajectories for viscosities η = 1.0 μPa•s (a, c) and η = 3.0 μPa•s (b, d) with respect to various polarizations (ζ) for a Brownian particle being optically manipulated in a trapping potential (a, b, corresponding to Fig. 1) or a saddle potential (c, d, corresponding to Fig. 3) are presented.The black dots indicate the coordinate origins and the starting positions of the trajectories, and the initial velocity of the particle is 0.